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Abstract 

In the study, the measurements of Scots pine height increments were used to compare the increments of pine trees 
of different age classes. All of the analyzed trees were growing in stands located on fresh mixed coniferous forest 
sites. The study concerned a 10-year period of growth of 8 tree age classes. Due to variation in climate conditions, all 
trees were studied over the same calendar period. Longitudinal analysis was used to compare different age classes 
of trees with reference to the increments in height. This procedure had not been previously used for such purpose. 
The results obtained did not confirm the hypothesis of parallel profiles implying that there existed differences in the 
growth of trees in various age groups. 
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Introduction 

Knowledge about relationships between the age of trees 
and their height increments is very important in man- 
agement of a forest system. Moreover, periodical height 
increments of trees are regarded as the major dendro- 
metric characteristic with direct influence on the ac- 
curacy of determination of tree and stand volume in- 
crements. The magnitude of tree height increments is 
needed when using methods to determine stand produc- 
tivity. The height increment is a variable characteristic 
and depends on several factors. The most important of 



these are tree age and biosocial position occupied by the 
tree, as well as soil and climatic or meteorological con- 
ditions (Beker 1997; Kazmierczak 2005; Zielski 1997). 

The regularity of the process of growing is inde- 
pendent from tree species (Assmann 1970). Differ- 
ences among species are connected with the extent and 
length of tree natural growth phases. In the first years 
of a tree life - the youth or juvenile phase - the incre- 
ment in height is not large. The next (second) period 
is called the height season or strength phase or vigor 
phase. During this time trees display the largest annual 
height increments, and the culmination point - the year 
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with the largest height increment - is observed. After 
a few years of large height increments, trees pass into 
the senescent phase: the period of stabilized increments, 
often much smaller than those produced in the second 
phase (or else with same size as before). 

The aim of this study was to analyze differentiation 
of annual height increments conditional on the age of 
trees. For this purpose, the statistical method of analy- 
sis of longitudinal data (also known under the name of 
profile analysis or longitudinal analysis) was used (mul- 
tivariate analysis for repeated measurement). Generally, 
the height increments in trees can be divided into the 
categories reliant on tree age. In the literature, there are 
available examples on the use of the data on trees in 
profile analysis. In the paper by Mehtatalo (2005), the 
longitudinal analysis was taken into account by esti- 
mating the models as random effects models using two 
nested levels: stand and measurement occasion level. 
Lukas and Diggle (1997) used the longitudinal analysis 
of repeated measurements to investigate relative differ- 
ences in the growth of Sitka spruce and Norway spruce 
seedlings during summer exposure to ozone over three 
growing seasons. In the study by Lappi (1997), the lon- 
gitudinal analysis was used in the simultaneous statisti- 
cal analysis of the height curves. 

Materials 

The experimental material included selected Scots 
pine stems from eight stands. The trees were of dif- 
ferent age (24, 33, 55, 62, 72, 84, and 92). The Draudt 
method was used as the selection method. In each age 
group, 25 trees were measured. Next, for each tree, 
ten measurements of increments were carried out over 
a ten-year growth period. This way, at 24 years of age, 
the first measurement actually applied to 15-year-old 
trees whereas at 33 years of age - the first measure- 
ment applied to 24 years of age, etc. All of the stands 
observed were located on fresh mixed coniferous for- 
est sites within the area of the Experimental Forest 
District Zielonka (temperate climate, E 17° 6’ 48.996”, 
N 52° 33’ 13.752 ). The same calendar growth period, 
extending from 1989 to 1999, was taken on for each 
studied tree, with the purpose to eliminate the effects 
of additional factors (e.g. climate and weather condi- 
tions) on the height increment. 



Methods 

The height was measured annually - always in the same 
trees, and thus the record values cannot be regarded as 
independent random variables, since the measurements 
were carried out on same experimental units. Longitu- 
dinal analysis was carried out with repeated measure- 
ments (performed on the same treatment at intervals) 
used to analyze growth variability. The applied model 
is: Y = XH + e, where Y (k x p) is the observation ma- 
trix. The columns of matrix Y are: y = [y^, ..., y jp \ 
for observation (trees) i = 1, 2, ..., n and group (age) 
j = 1, 2, ..., k , where p is the number of time points 
(/ = 1, 2, . . ., p). Then X is the ( nk x k ) design matrix, S 
is the (k x p) matrix of unknown parameters and e is the 
(nk x p ) error matrix built from the measurement errors 
e jX . In the above model it is assumed that e/ ; - = \e ijX , ..., 
e //p \ are independent, and have the multivariate normal 
distribution with expectation equal zero and unknown 
nonsingular, full column rank, covariance matrix 2, 
which is identical for each i = 1 , 2, ..., n,j = 1 , 2, ..., k. 

Three hypotheses were tested in this analysis. The 
first was the hypothesis H ol which assumed the ab- 
sence of interactions between calendar years and stand 
age, i.e. parallelism of the profile of height growth for 
trees in different age groups. If there is no basis to 
reject the null hypothesis (growths of height in the 
studied period were independent of tree age), then 
the hypothesis H 02 assuming equality of the growths 
(profiles) is checked, as well as the hypothesis H m as- 
serting the absence of changes in the time of medium 
growth for stands (in other words: absence of a trend, 
equal growths in every year of the studied period, ab- 
sence of differences between years). For verification of 
the first and the third hypothesis, Roy’s largest root is 
generally used, while the second hypothesis is tested 
using the Hotelling T 2 statistic (Morison 1990). If the 
hypothesis H m is rejected, Roy’s (or Scheff’s) simul- 
taneous confidence intervals can be used to find the 
age groups of trees which contributed to this. Equality 
of the expectation for different study years should be 
tested in the same way. In this case T 2 statistics for de- 
pendent variables are used separately for each group. 
If the hypothesis H ol is not rejected but H 02 and/or H m 
are, then Scheff’s intervals can be used to find the rel- 
evant differences between groups and time points. The 
details of the method used were described by Zawieja 
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(2012), however, in the present study, the increments 
of trees in the same age group were checked among 
different Kraft classes. All the conclusions were based 
at a significance level p = 0.05. 

Results 

Before attempting the analysis, it is worth to present 
general trends (profiles) of groups’ means with respect 
to the time points. In figure 1, there are presented the 
profiles for the tree age groups studied. The increments 
of the youngest trees are most changing in the years, 
while the increments of the oldest tress are most stable. 
It should be noted, however, that increments of all age 
groups have a similar pattern, but the youngest tree, 24 
and 33, strongly react to weather conditions (for exam- 
ple 1993, when precipitation was three times higher than 
in other years - the increments were smaller). Thus, it is 
reasonable to examine the hypothesis concerning paral- 
lelism of the profiles of individual groups. 

Based on profile analysis, it was found that the 
hypothesis H m (profile parallelism) ought to be reject- 
ed, the Roy’s statistic: 9 = 2.10997, the critical value: 
6 mb = 0.152 with parameters S = 7, M = 0.5, N = 80.5 
(Zawieja 2012). Rejection of the aforementioned hy- 
pothesis entailed rejection of the hypothesis of profile 
equality as well as that of equality of time -points ef- 
fects. 

Because the hypothesis of parallelism of profiles 
was rejected, it was necessary to check which age 



groups contributed to this non-parallelism of profiles. 
On the basis of the Scheff’s confidence intervals it was 
concluded that the height increments of the youngest 
(24) and the oldest (92) trees were clearly different from 
other groups. The trees in age groups from 55 to 84 
formed the largest homogeneous group (tab. 1). 



Table 1 . Multivariate comparison of multiple groups of trees 
(stands) at significance level p = 0.05 



Age group 


24 


33 


43 


55 


62 


72 


84 


92 


Homogeneous 

groups 




a 


a 
















b 


b 
















c 


c 


c 


c 





The hypothesis of equality of time-points H 0 3 was 
rejected, thus calendar years were compared to indicate 
the year or years in which the increments of all groups 
differed significantly. In this case the pairs of consecu- 
tive years (the first year with the second, second with 
the third, etc.) were compared. No significant differ- 
ences between 1995 and 1996 and also between 1996 
and 1997 were shown (tab. 2). 

Since the hypothesis of the equality of profiles was 
rejected, ten one-dimensional analyzes of variance were 
performed separately for each study year. The following 
values of F were obtained (in subsequent groups, from 
the youngest to oldest) 100.58, 78.84, 202.72, 181.56, 
59.80, 117.29, 81.36, 49.66, 46.45, 40.62 with seven de- 
grees of freedom for the nominator and from 164 to 167 
(because the outliers were removed) for the denomina- 






Figure 1 . Mean Scots pine height increments in observation years and profiles of age-group trees in ten-year period 
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tor, /7-values were always less than 0.0001. In all of the 
analyzes, the hypotheses of equality of tree age groups 
was rejected. Hence, the simultaneous multiple com- 
parisons were performed. The results of the compari- 
sons are shown in table 3 and it should be noted that 
the height increments in the group of the youngest trees 
differed from the increments in other age groups (from 
the age 15 to 22). In 1997, when the trees reached the 
age of 23 years, they did not differ from the trees from 
the age groups 43 and 33. The remaining age groups of 
trees were combined into larger or smaller homogene- 
ous groups with reference to the height increments. The 
trees from the age groups: 55 and 43 years old, moreo- 
ver - those from 62, 72 and 84 years and also 72, 84 and 
92 years formed a homogeneous group in almost every 
study year (with the sole exception in 1991 for the age 
groups: 62, 72 and 84 years). 



Table 2. Comparisons (contrasts) between successive 
calendar years of measurement (time points) 



Study year 


F-test 


Significant difference 


(time point) 


(p-value) 


1989-1990 


121.4 


<0.0001 


1990-1991 


158.3 


<0.0001 


1991-1992 


56.6 


<0.0001 


1992-1993 


236.6 


<0.0001 


1993-1994 


107.12 


<0.0001 


1994-1995 


162.0 


<0.0001 


1995-1996 


0.21 


0.6512 


1996-1997 


0.06 


0.8040 


1997-1998 


4.76 


0.0305 



Table 4. Comparisons (contrasts) of successive years of 
measurement in individual stands (tree age groups) 



Age-group 


24 


33 


43 


55 


72 


84 


1989-1990 


- 






- 


- 


- 


1990-1991 






- 




- 


- 


1991-1992 




- 


- 


- 


- 


- 


1992-1993 


- 








- 


- 


1993-1994 


- 




- 


- 


- 


- 


1994-1995 


- 


- 


- 









**** significant differences at p = 0.00001; - no significant 
differences. 



The tree age groups were highly different from each 
other as shown above. From here, separately for each of 
the age group, profile analyses were performed. Con- 
sequently the value of T 2 statistics exceeded the criti- 
cal values (value of T 2 statistics for subsequent groups, 
from the youngest to oldest, were following as follows: 
460.17, 263.68, 318.45, 315.67, 160.53, 340.98, 176.68, 
207.31, the critical value of Tl 05 . gA5 = 5.853; all ^-values 
were < 0.0001) for every stand, so in each group the null 
hypothesis of absence of differences between the years 
was rejected. Scheff’s confidence intervals were used for 
comparisons between increments in the successive study 
years (for each contrast). The results of the comparisons 
are given in table 4 (only for these age groups and pairs 
of years when significant differences existed). In the age 
groups: 62 and 92 years, significant differences between 
successive years was not detected. The height increments 
frequently differed between the years 1990 and 1991 (age 
group: 24, 33, 55 years), similarly in 1994 and 1995 (age 
group: 55, 72, 84 years). For the respective age groups, 



Table 3. Homogeneous groups of stands (tree age determined in the year of measurement) in successive calendar years of 
measurement (p = 0.05) 



1989 


1990 


1991 


1992 


1993 


1994 


1995 


1996 


1997 


1998 


15 


16 


17 


18 


19 


20 


21 


22 


23 a 


24 a 


24 a 


25 a 


26 


27 


38 a 


29 


30 a 


31 a 


42 ab 


43 ab 


34 ab 


47 ab 


48 a 


49 a 


28 a 


39 a 


40 a 


41 a 


32 ab 


33 ab 


46 b 


35 b 


36 a 


37 a 


50 a 


51 ab 


52 a 


53 ab 


54 be 


55 be 


53 c 


54 be 


55 


56 b 


57 ab 


58 be 


59 b 


60 be 


61 cd 


62 cd 


75 cd 


64 cd 


65 b 


78 be 


67 b 


68 cd 


69 be 


70 bed 


71 de 


84 cd 


63 cd 


76 cd 


77 b 


66 be 


79 b 


80 cd 


81 be 


82 cd 


83 de 


72 d 


83 d 


84 d 


85 b 


86 c 


87 b 


88 d 


89 c 


90 d 


91 e 


92 d 
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there were found clear differences in the average incre- 
ment size between the years, which is presented in figure 
1. Similarly, significant differences in the increments of 
the age group 33, between the pairs of years: 1989 and 
1990, 1990 and 1991, 1992 and 1993 as well as 1993 and 
1994 can be seen in figure 1 showing average increments. 

Discussion 

In the present study, profile analysis method was used 
to determine differentiation of the height increments 
depending on Scots pine age classes. This is a new ap- 
proach to analyzing this problem. For example, in the pa- 
per by Crecente-Campoa et al. (2010), an individual-tree 
growth model was developed. Stankova and Dieguez- 
Aranda (2013) presented modeling height increments 
based on other characteristics of trees. Gerelbaatar and 
Baatarbileg (2011) studied height increment changes in 
close correlation with diameter increment for different 
vegetation zones. Profile analysis do not give a model of 
growth but on the basis of this analysis it can be seen 
if the increments in the different groups of trees are the 
same or not in the successive years. If the increments are 
not the same, it can be concluded which groups differ 
from each other. As mentioned above, the results of pro- 
file analysis provide information as that obtained in the 
analysis of growth. Namely, one can conclude at which 
age the increments of trees do not change from year to 
year, in the analysis of growth, there is often observed 
autocorrelation, in profile analysis the effects of autocor- 
relation and weather conditions are excluded, in view of 
the fact that the analysis is tailored to data concerning 
repeated measurement on same experimental units. The 
used method of analysis for repeated measures seems to 
be a very good method for evaluating the increments of 
tree heights, and what is more, it can be used for other 
than Scots pine tree species as well. 

Conclusions 

Calculations were performed by means of profile analy- 
sis allowed to draw the following conclusions: 

- The height increment of trees representing stands of 
the age groups 24 and 92 is different from the height 
increment of trees from other age groups. 



- Additionally, it was clear that the youngest trees 
(24-year-old) finished the vigor phase of their height 
increment at the age of 22 years and then they passed 
into a less intensive growth stage. 

- The results obtained indicate that in the year 1991, 
the homogeneous groups in terms of tree height in- 
crements differed from those observed in adjoining 
years. 

- in the year 1991 there are only two homogenous 
groups of younger and older trees, missing interme- 
diate groups linking these two groups (see tab. 3). 
This suggests that the conditions in that year were 
different than in adjoining years. Because the cause 
of it is not visible (weather conditions do not dif- 
fer noticeably from other years, as it visible in 1993 
- heavy precipitations) thus this problem will be in- 
vestigated in further research. 
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